Electron spin resonance in high-field critical phase of gapped spin chains 



A. K. Kolezhukfcl 

Institute of Magnetism, National Academy of Sciences and Ministry of Education of Ukraine 
36(B) Vernadskii avenue, 03142 Kiev, Ukraine 
Institut fur Theoretische Physik, Universitat Hannover, Appelstrafie 2, 30167 Hannover, Germany 



o 
o 

(N 



H.-J. Mikeska 

Institut fur Theoretische Physik, Universitat Hannover, Appelstrafie 2, 30167 Hannover, Germany 

(February 1, 2008) 

Motivated by recent experiments on Ni(C2HgN2)2Ni(CN)4 (commonly known as NENC), we 
study the electron spin resonance in the critical high-field phase of the antiferromagnetic S = 1 
chain with strong planar anisotropy and show that the ESR spectra exhibit several peculiarities in 
the critical phase. Possible relevance of those results for other gapped spin systems is discussed. 
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I. INTRODUCTION 



Recently, there has been a growing interest in the 
properties of low-dimensional r^ipi systems subject to 
strong external magnetic field. EfO In gapped spin sys- 
tems, when the external magnetic field H exceeds a crit- 
ical value H c , it closes the gap and drives the system into 
a new critical phase with finite magnetization and gap- 
less excitations. When the field is further increased, the 
system may stay in this critical phase up to the satura- 
tion field H s , above which the system is in a saturated 
ferromagnetic state. Under certain conditions, however, 
the excitations in this high-field phase may again ac- 
quire a gap, making the magnetization "locked" in some 
field range; this phenomenon is known as magnetiza- 
tion .plateaus and has been receiving much attention as 
welLBta 

Recently, electron spin resonance (ESR) experiments 
on Ni(C2HgN2)2Ni(CN)4 (commonly abbreviated as 
NENC) in strong magnetic fields were conducted.LJ This 
compound is believed to be a realization of the S =p4 
chain with strong planar and weak in-plane anisotropy.E-D 
The theory describing ESR response for this system out- 
side the critical phase (i.e., H <. H c or H > H s ) was 
developed by Papanicolaou et al.0 However, much of the 
ESR data of Ref. ^ belongs presumably to the field range 
H c < H < H s , i.e. exactly to the region where the the- 
ory is lacking, which makes the interpretation of the data 
rather difficult. 

Motivated by those experiments, in the present pa- 
per we study theoretically the zero-temperature ESR re- 
sponse in the critical phase of a planar S = 1 chain. 
It is shown that the typical feature of ESR in the criti- 
cal phase is the appearance of continua with resonances 
being determined by power-law singularities instead of 
well-defined quasiparticle peaks. We predict that a char- 
acteristic change of the slope in the field dependency of 
the ESR resonance frequency takes place at the critical 
field H — H r . 



The paper is organized in the following way: in Sect. 
H we introduce the effective Hamiltonian for the planar 
chain, Sect. [II is devoted to the calculation of resonance 



frequencies and exponents characterizing the correspond- 
ing singularities, and, finally, in Sect. IV we discuss pos- 



sible existence of similar features in other gapped one- 
dimensional spin systems. 



II. PLANAR S = 1 CHAIN: EFFECTIVE MODEL 

We start from the model of strongly anisotropic anti- 
ferromagnetic 5 = 1 chain described by the Hamiltonian 



H 



?z\2 



n n 



(i) 



where the planar anisotropy D is assumed to be much 
stronger than the exchange constant J, E <C D is a weak 
in-plane anisotropy, and h = gfJ-sH , where g is the Lande 
factor, \ib is the Bohr magneton, and H is the external 
magnetic field. The in-plane anisotropy should be in- 
cluded since it spoils S z as good quantum number and 
thus allows certain transitions which are otherwise forbid- 
den; generally, the in-plane anisotropy constant E may 
be comparable to J. 

Let us consider first the noninteracting case J = 0, 
then (Q) amounts to a single-ion problem. The spectrum 
of a single ion consists of three states, whose wave func- 
tions and energies read as follows: 

|») = |0), e v =0, 

| a) = cosa|+) — sin a | — ), e a = D — h 

\b) = sin a\+) + cos ct| — ), e& = D + h 



h=(E 2 + h 2 ) 1 ' 2 , 



sina = —=(\ - h/h) 1/2 , 
v2 



(2) 
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where |0), |±) denote the spin-1 states. For J <§; 13 and 
weak field h < h c , the ground state of the model (|l|) 
can be described, to a good approximation, as a direct 
product of \v) states of separate magnetic ions. The el- 
ementary excitations are propagating \a) and \b) states, 
the spectrum has a gap, and the lowest excitation is a 
degenerate doublet. This doublet gets split off by the ex- 
ternal field, | a) state coming down with the field, and \b) 
states going up. The low-temperature ESR response for 
h < h c is determined by the transitions fepm the ground 
state of the type \v) — ► \a) and \v) — ► \b)E 
When the field h exceeds the critical value 

h c = {{D -2J) 2 - E 2 } 1/2 , 

the system enters critical phase with a finite density M 
of | a) states, and a new type of ground state transitions, 
namely of the |o) — * \b) type, becomes possible (note 
that this transition is allowed only in presence of the fi- 
nite in-plane anisotropy E). The "old" types \v) <-> \a) 
and \v) — > \b) still remain possible. As the density M 
increases from to 1, the \v) — ► \b) signal should become 
weaker, while the intensity of the \a) — > \b) transition 
should increase, 
mined by the equation 



$afc y^(a w 6 n +i + b n a n+ i + h.c), 



Above the saturation field h s , deter- 



h = D + 2j{l + (h/h) 2 }, 

the density of a-particles is equal to one. 

In order to describe this picture quantitatively, we first 
introduce the hardcore boson operators (V n , b' n creating 
respectively the |o) and \b) states from the vacuum state 
\v). Not more than one boson is allowed to be present 
on any site, which defines the set of physical states. The 
effective Hamiltonian in terms of hardcore bosons can be 
written as H c s — V(Ho + 'Hi U _t)V, where V is the projec- 
tor onto the set of physical states. The "unperturbed" 
Hamiltonian 7io has the following form: 

Ho = 53(A«4an + A 6 6j l 6 n ) 

n 

+ 5Z^ Q a n a n+l + *b b l b n+l + h.C.), (3) 

n 

where the self-energies are A a ,h — e a .b and the hopping 
amplitudes t a = tb = t = J . The interaction part Hint 
looks as follows: 

Hint = U aa '^2(al l a n )(al l+1 a n+1 ) (4) 

n 

+ U ab J2{( a l a n)(Ki+i b n+i) + (bib n ){al l+1 a n+ i)}, 

n 

where the interaction constants are given by U aa — 
-U ab = J(h/h) 2 . 

Our goal is to take into account the effects arising in 
the first order in the coupling J, and therefore we have 
neglected the pair creation term 



with $ a h = J (h/h), whose contribution to any energy 
level arises only in the second order in J. This simplifies 
the problem considerably, making the total numbers _/V a 
and N b of a and b particles good quantum numbers. At 
the same level of approximation, one can neglect the pro- 
cesses of exchanging the positions of neighboring a and b 
particles, described by the term 

t ab ^2 {blal +1 a n b n+1 +h.c.}, 



since the 'interparticle' hopping amplitude t a b — 
J(E/h) 2 is very small in the critical region, E/h ~ 
E/D <C 1. Further, we will be interested in configu- 
rations with at most one b boson, so that the interaction 
between b particles can be also safely neglected. 

Let us first consider the "unperturbed" Hamiltonian 
VHq'P (note that the physics described by this Hamil- 
tonian is nevertheless nontrivial due to the single occu- 
pancy constraint). 

In absence of 6-particles (Nb — 0) the spectrum of 
the problem can be in principle obtained by mapping to 
noninteracting fermions with the help of the well-known 
Jordan- Wigner transformation, but it is more convenient 
to stick to the hardcore boson language. From the Bethe- 
ansatz solutionlla one knows that the energy spectrum is 
given by 



E({k t }) = J2{A a + cos kt), 



2tt , s 

+ T h, (5) 



where the numbers Ii are all different and respectively 
half-integer (integer) if the total number of particles 
N = N a is even (odd), and L is the number of sites 
which will be assumed even for convenience. The total 
momentum P of a state defined by a certain choice of 
N different numbers /; is P = J2i h- The ground state 
|g.s.) is given by a symmetrical dense distribution of /; 
around zero: 



h 



N-l N-'S 



N -1 



(6) 



describing the Fermi sea of particles with momenta in the 
interval [hp, 2tt — hp], where we shall define 



7r(l -N/L). 



(7) 



The total momentum of the ground state configuration 
is P = (mod 2ir) if N is even, and P — ir if N is odd. 

Since the hopping amplitudes for a- and 6-particles are 
equal, it is easy to realize that the above picture of the 
distribution of wave vectors remains true when one has 
the total number of particles N = N a + Nb, N a of them 
being of the a type and Nb of the b type: they form a 
single "large" Fermi sea. 
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III. ESR TRANSITIONS 



We would like to study the ESR transitions from the 
ground state, which survive in the low-temperature limit. 
The ESR intensity I(ui) for the ground state transitions 
is determined by the formulalij 



J(w) oc a; X) I EW^StlS' 8 - 



5{E f ~ E g . s . - u), 



(8) 



where h" are the components of the radio frequency field, 
and / labels all possible final states. 



where Ai is the normalization factor, and the phase 
shift Si = S = —ir/2 is independent of the state la- 
bel I in our case of noninteracting hardcore particles. 
The factor (— ) s = ±1 is the sign of the permutation 

l> v> ' '.'.$)> where & < 4 < ••• < y'i N > ^ich 

ensures that the total wave function is symmetric under 
the permutation of any two particles. The allowed values 
for the wave vectors k[ of the excited state are given by 
(ft) with integer // . The energy of this state is given by 



i 



(A a + 2t cos k'i) + A fc + 2t cos A, 



(11) 



A. \v) i— » \b) transitions 

Let us first consider the \v) i— ► \b) process with AN a = 
and ANi, = 1, as the most interesting one (note that it 
is allowed even in absence of the in-plane anisotropy E) . 
This process is determined by the operator S 1 ^, which 



is in this case equivalent to y 1 — h/hJ2 n Hi- Let us as ~ 
sume for convenience that the initial number of particles 
N = N a is even, so that N = 2N . Changing the total 
number of particles by one causes, according to (Q), the 
rearrangement of the allowed wave vectors. This leads 
to the vanishingly small transition matrix elements in 
(||) (orthogonality catastropheEJ) , however it is compen- 
sated by a diverging number of states |/) having nearly 
the same energy. The problem of calculating the ESR 
response in this case amounts to that of calculating the 
spectral properties of a mobile hardcore impurity sud- 
denly created in a hardcore boson system, the impurity 
having the same mass as the other particles, and the dis- 
persion law being of the type A a & + 2tcosfc. A similar 
problem was studied in Refs. |2l|-|24] and is closely related 
to the so-called Fermi edge singularities in the photoemis- 
sion/adsorption spectra (see, e.g., Ref. and references 
therein). 

The excited wave function |/) can be generally repre- 
sented in the following form: 



l/) = 



$(y'u--,y'N)\ x 'o; vi-'-vn), (9) 



where x' denotes the position of the b particle and 
x\ = x' Q + y[, i = l,...,N denote the positions of the 
a-particles (so that y[ are their coordinates relative to 
the 6-particle) , and P' is the total momentum of the ex- 
cited state. !_. 

Following Castella and Zotosj£3 one can write the re- 
duced wave function $ in the dctcrminantal form: 



det{v>i (?/■)}, 



A' 



<p,(y) = A l [e^y+ S ^ -jjY, e i{k '~ V+5n) ), (10) 



where A is defined by 



(12) 



and plays the role of the momentum of the b particle. 

It is easy to see that the wave function 6|. Jg.s.), ob- 
tained from the ground state by creating a 6-particle at 
site Xq, can be written in a similar determinantal form: 

&t o | g . s .) = e iPx ° Y, f(Vu.--,VN)\xo\Vi...VN), (13) 



where yi = x% — xq are again relative coordinates of a- 
particles, P = ki is the total momentum, and the 
reduced wave function 



f(yi, ■ ■ ■ ,vn) = (-) J 



i 



det{i/ji{yj)}, 



(14) 



is built from simple plane waves ipi(y) — (1 / ' \[V)e %klV , 



yi V2 



{— ) s is the sign of the permutation ( 

where y h < y l2 < ■ ■ ■ < y iN . 

Now, one can see that the matrix element entering (g) 
takes the form 



(f\YbU&*-)=^LS(P,P')M fi , 



<Vn) /(2/i,2/2,. ..,2/iv), (15) 



V1V2 — VN 



where 6(P, P') is the Kronecker symbol telling us the 
well-known fact that the total momentum should not 
change during the ESR transition, P = P' (mod 2ir). The 
reduced matrix element Mfi can be represented as a de- 
terminant of the overlap matrix between 'old' and 'new' 
wave functions, 



M fi =det{(HVz)}. 



(16) 



As shown in Ref. ^2|, the wave functions fi{y) defined 
in (fil]) behave asymptotically as plane waves, fi{y) i— > 



3 



<pi{y) = (l/VL)e i( - k ' y+5 '\ and the overlaps (||) calcu- 
lated with ipi can be with the accuracy of 0(hxL/L) re- 
placed by those calculated with fy. Then the matrix 
element Mfi is asymptotically equal to the overlap be- 
tween two Slater- type wave functions, one describing the 
system with the total momentum P = Jbj an d the 
other corresponding to the system with the momentum 
Q = J2i k'i = P 1 — A. It is thus clear that Mfi can be 
nonzero only if P = Q. This gives us the complete set of 
selection rules as 



P = P'. A = 0. 



(17) 



The ground state configuration is given by the follow- 
ing distribution of the wave vectors: 

fc, = 7r ±^(l + ~), l = 0,...,(N -l) 

and has the total momentum P — 2ttN q , while in the 
excited state the numbers Jj are integer, and the lowest 
energy configuration built according to (^) would have 
the total momentum / k(2Nq + 1). In order to get the mo- 
mentum change AP back to zero, one has to introduce 
some excitation. In the simplest way this can be achieved 
by creating a hole at k = it, i.e., one gets the symmetric 
excited state configuration |/ s ) given by 



2tt . 
7r ± — (n - 
L 



0. 



,N -1. 



(18) 



The overlap Mj { can be easily calculated by means of 
passing to the even/odd states (i.e., to sines and cosines 
of kiy and (My + S), respectively). The determinant 
then factorizes into a product of two (in our case equal) 
Cauchy-type determinants M", 



M (±) = dct 



1 



(n — n') — S . 



Determinants of this type were calculated by AndersonEj 
and shown to be algebraically vanishing in the ther- 
modynamic limit, M^' cx L~P ± / 2 , with the exponents 
f3± = (S/tt) 2 = 1/4, so that 



\M s fl \ 2 cx L-A 



(19) 



with the orthogonality catastrophe (OC) exponent f3 s — 
j3+ + /9_ = 1/2. This exponent can be also calculated in 
a different way, using-the results of boundary conformal 
field theory (BCFT).c3 For this purpose it is necessary to 
calculate the energy difference AEf between the ground 
state and the excited state |/), including the 1/L correc- 
tions. Then in case of open boundary conditions the OC 
exponent /3, according to BCFT, can be obtained as 



ILAEf _ 2AE f 



TTVp 



AE m \r 



(20) 



Here vf = 2tsinfcF is the Fermi velocity, so that 
A-E m i n = Tivp/L is the lowest possible excitation energy, 



and AEf is the 0(1/ L) part of AEf (i.e., with the bulk 
contribution subtracted). In this last form this formula 
should be also valid for the periodic boundary conditions, 
then Ai^min should be replaced by 2itvf/L. It is easy to 
calculate AEf for the symmetric configuration |/ s ): one 
gets AE s f = A h + 2t + AE' S , where 



N -l 



. / f 27T,, 1, 27T,, ,1 

= At { cos —(/ + -) - cos — (I + 1)} 

1=0 

= 2t(l + cos k F ) + ty smk F + 0(1/ L 2 ). (21) 

Applying (^) , one gets the correct value = 1/2 which 
coincides with that obtained through a traditional way 
by calculating determinants. 

The symmetric excited state configuration with Q = P 
is obviously not the configuration with the lowest energy. 
The symmetric configuration is, however, the configu- 
ration with the highest overlap with the ground state. 
Indeed, there are configurations with asymmetric distri- 
bution of momenta around k — n, which nevertheless 
satisfy the selection rules, e.g., the following one which 
we will denote as |/ a (n)): 

2ir 

k[ = it + —l, l = -(N -n),...,(N +n-l). (22) 

If we require that Q = P + 2t:A, where A is an integer, 
the following condition on hp is obtained: 



2m- 1 
2n - 1 ' 



A. 



(23) 



which gives the dense set of allowed values of hp, at which 
the configuration of the type p2] ) will satisfy the selec- 
tion rules. The energy difference from the ground state 
energy is AE a f (n) = A b + 2t + AE' a (n), with 



AE' a {n) =t(2n-l) 2 -sinfc / + (9(l/ J L 2 ) 



L 



(24) 



The corresponding OC exponent (3 a (n), according to 
©; is given by 



A» = ^(2n-1) 2 



(25) 



and is always larger than the OC exponent for the sym- 
metric configuration j3 s — 1/2 (except for n = 1, which 
would mean kp = it, i.e., the vanishing Fermi sea). 

The main contribution into the ESR response is thus 
coming from the states whose energy is close to that 
of the symmetric configuration. The summation over 
"shake-up" (i.e., particle- hole) excitations around the 
symmetric configuration leads Jm the singularity in the 
ESR intensity I(oj) of the formEa 



I(u>) oc l/(u> — Wq)", uj > LOq 



(26) 
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with a = a s = 1 — /3 S = 1/2, and the threshold fre- 
quency wo = AEj = A b + At + 2t cos kp. It is remarkable 
that the contribution from the asymmetric configuration 
displays no singularity, since the corresponding OC ex- 
ponent j3 a > 1, so that a a < 0; this contribution yields 
just some scattered intensity at frequencies u> > A 5 + 2t, 
see Fig. [|. The integrated intensity of \v) — > |a) contri- 
bution is proportional to the total number L(l — M) of 
empty sites, and thus decreases with the magnetic field; 
this contribution does not exist above the saturation field 
H s . 

Taking into account the definition of kp, namely A a + 
2t cos kp = 0, the formula for threshold ojq can be rewrit- 
ten as 



A a (M) + 2t cos k F = 0. 



(29) 



uj (v -> b) 



At 



(27) 



One may notice the following remarkable feature: for 
weak magnetic fields below h c the \v) 1— » |6) process leads 
to a quasiparticle peak in the ESR response at the fre- 
quency u>' Q — Ah + 2i, so that below the critical field the 
slope du}' /dh of the ESR line as a function of the field 
h is approximately 1 (for E <C D). For h > h c the 
same process yields the Fermi-edge-type singularity ( p^ ) 
in the ESR response, with the edge frequency loq being 
given by (|27]). If one makes a reasonable assumption, 
that experimentally observed ESR line for h > h c fol- 
lows the behavior of the edge singularity, we get for the 
slope of the ESR line above the critical field the value 
dcoo/dh w 2. Thus, the slope of the experimentally ob- 
served line should change abruptly at the critical field 
h = h c , as shown schematically in Fig. |l|. It is worth- 
while to note that such a behavior is .reminiscent of the 
picture experimentally seen in NENCQ in the intermedi- 
ate field regime. 

Though taking into account the effects of temperature 
is beyond the scope of the present paper, we would like 
to remark that the singular contribution of \v) — > \b) pro- 
cesses to the ESR response should be enhanced at finite 
temperature, due to the presence of a finite number of 
holes in the ground state. The density of states has a 
singularity at the bottom of the Fermi sea, so that the 
main contribution will come at the same frequency (27). 



This could be important for interpreting the results of 
experiments on NENC,Q where the temperature was of 
the order of J. 

If the interaction terms (^) are taken into account, the 
corresponding effective model cannot be solved exactly 
any more, and we cannot calculate the OC exponents in 
this case. However, one can treat the effect of interac- 
tions in a sort of the "mean- field" approximation, i.e., 
replacing simply 



A a h+ A a {M) = A a + U aa M, 
A b ^ A b (M) = A b + U ab M, 



(28) 



where the particle density M = A — kp /tt has to be cal- 
culated self-consistently from the equation 



One can check that such an approximation, though being 
crude, delivers correct values for the critical fields H c and 
H s . Actually, in absence of o-particles one can show that 
the effect of U aa is not only to change the self-energy A a , 
but also to renormalize the hopping amplitude t a (in the 
first order in U aa one gets t a 1— » t a + (U aa /n) sinfcp), and 
both t a and U aa are of the order of J, so this correction, 
is not negligible. According to the numerical studies,EZl 
the OC exponent (3 is rather sensitive to the ratio of 
hopping amplitudes r — t b /t a (the behavior of (3(r) is 
approximately linear around r — 1), so that the corre- 
sponding OC exponents will certainly change compared 
to the noninteracting (hardcore) case, and they will be- 
come fcp-dependent. One may nevertheless expect that 
at the qualitative level our picture of transitions remains 
true, particularly, the effect of changing the slope of ESR 
lines at H = H r should survive. 



B. \a) 1 — ► \b) transitions 



The I a) 1— > |6) processes are possible only in presence 
of the finite in-plane anisotropy E, which mixes |— ) and 
|+) states into \a) and |6), respectively. Those transitions 
are determined by the operator S^ ot , which is in this case 
equivalent to {E/h)Yl n b\i a n- One can easily see that 
now the wave functions Y^ n 6iOjj|g.s) and |/) can be both 
represented in the form (^), ([L0|) with N — > [N — 1). The 
total momentum of the ground state can be represented 
as P — X^iLi 1 ki + \ an d the momentum of the excited 
state has a similar form, P' = X^i 1 + Through the 
same line of arguments as before one gets the selection 
rules 



P = P', A = A. 

Since the number of particles was not changed during the 
transition, the allowed values of wave vectors k[ are now 
the same as fcj and are given by 7r + (2tt /L)Ii with half- 
integer Ii. Thus, there is no orthogonality catastrophe in 
this case, and one gets a well-defined quasiparticle peak 
at the frequency 



A„ 



(30) 



One can visualize this process as simply replacing one of 
the a-particles in the Fermi sea of the ground state by 
a 6-particle, without changing the wave vector distribu- 
tion. The corresponding ESR line goes parallel to the 
\v) — > 16) line, slightly below it, as shown schematically 
in Fig. |l]. The intensity of this line is proportional to the 
total number LM of a-particles and increases with the 
field; this line continues to exist at H > H*. 
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C. \v) <-> | a) transitions 

Those transitions are allowed in absence of the in- 
plane anisotropy. The number of particles changes by 
one, which leads to a change in the distribution of mo- 
menta and indicates absence of the quasiparticle peak. 
Since only two states per site are involved (iVf, remains 
zero), those processes can be treated within the effective 
spin-i model, similarly to the way it was done for spin 
ladders Bli-3 If one neglects the small corrections of the or- 
der of (E/h) 2 , this effective model is the easy-plane XXZ 
chain in an effective longitudinal magnetic field with the 
HamiltonianEj 



I(ui) oc 



H, 



eff 



E 



BS* 



where 5 M are the spin operators of the effective spin-i 
chain, J xy — 2 J, the easy-plane anisotropy J z j J xy = 1/2, 
and the field B = h- J - D. 

Calculating the low-frequency ESR response amounts 
then to the knowledge of the dynamical structure factor 
(q, uS) of the effective XXZ chain, I(uS) oc (q = 
0,oj). The asymptotic behavior of the correlation func- 
tions fcpitke XXZ chain is known from the bosonization 
resultsHy 



(S+(x,t)S- (0,0)) 
C 2 



Ci 



• cos(7ra;) (31) 



(x 2 - v 2 F t 2 y/^ K ) 

gi(2/cp — tt)x g — i{2k,F—7r)x ■ 



-A 



(x 2 - v 2 F t 2 ) K + l /^K)-i I (x - v F t) 2 (x + v F t) 2 I ' 

Here C\ t 2 are some constants, vf is the Fermi velocity 
and K is the so-called Luttinger liquid parameter. For 
B = (which corresponds to kp = tt/2, or M = 1/2, or 
equivalently to the physical field H in the middle of the 
critical region, H = \{B C + H s )) and J z /J xy — 1/2 one 
has0 



v F 



= 3^3/2, # = 3/4. 



(32) 



At nonzero B (i.e., for M ^ 1/2) the values of vp and K 
can be obtained numerically by solving the Bethe ansatz 
equationspj the corresponding solutions are presented in 
Fig. 2 of Ref. [l2| For kp — > 7r (0), which corresponds 
to H — > H c {H s ), respectively, i.e., when the number of 
particles (holes) becomes small, K tends to 1/2 and vf 
to zero. 

Fourier-transforming (|3l|), one gets the q = dynamic 
structure factor (q = 0,u>) having edge-type singu- 
larities at the frequencies 



Wo(Q) = Vfl 



(33) 



where Q = it and Q = 2k p — 7r for the contributions of 
the first and second term in (|3l|), respectively. The sin- 
gularities are of power-law type, with the exponent being 
determined by the Luttinger liquid parameter K: 



•wo(Q)] 



v(Q) 



V {2kt 



K ~IK- ( ^ 



Thus, two lines corresponding to the singularities in the 
\v) «-> | a) processes at ui = ujq(tt) and u> = uo(2kF — 
7r) should be visible in the ESR spectrum, as shown 
schematically in Fig. |l|. Inclusion of the finite temper- 
ature would lead to damping of the singularities 



IV. DISCUSSION 

In the model of planar S — 1 chain considered 
in the previous section we have seen several features 
of the ESR spectrum in the high-field critical phase: 
appearance of low-lying and high-lying continua with 
power-law singularities, counter- intuitive change of the 
slope of ESR lines at the critical field H = H c , etc. 
One may ask oneself whether similar features can be 
present also in other one-dimensional spin systems, e.g., 
strongly coupled S = | ladders like (^^(CsH^N^^CLi 
(usually abbreviatecL-CuHpClH) and (CsH^^^CuB^ 
(abbreviated BPCBEl), or S = 1 Haldane chain 
Ni(C5H 14 N 2 ) 2 N 3 (PF 6 ) (known as NDMAPcffl) . 

As far as the model of S = i ladder with strong rung 
interaction Jr and weak leg coupling Jl <C Jr is con- 
cerned, one can construct the effective model exactly in 
the same way as for the large- 1? S = 1 chain. The role 
of vacuum state \v) will be played by singlet state |s) on 
a single rung, and one will have three types of "parti- 
cles" corresponding to the three triplet states \t±), \to). 
Applying a strong magnetic field will lead to closing the 
gap and to a finite density of \t + ) states in the ground 
state in the critical phase. However, in absence of any 
additional couplings (e.g., Dzyaloshinskii-Moriya interac- 
tion or non-uniaxial anisotropies) the only possible direct 
transitions from the ground state are determined by the 
processes of the \t + ) — > \to) type, which, similarly to the 
| a) — > | b) processes considered in the previous section, 
yield a quasiparticle peak with the resonance frequency 
cij + o = h. Only in presence of such additional interactions 
providing finite admixture of triplets in \v) one could see 
\v) ~ * l^o and \v) — * lines, which should then exhibit 
the change of slope at H = H c with the simultaneous 
appearance of the continuum above them. Similar argu- 
ments apply also to the S = 1 Haldane chains: generally, 
in order to observe the "interesting" lines, one needs to 
have some perturbations allowing the direct transitions 
from the ground state at H = 0. 

In this sense, the planar S = 1 chain is a remark- 
able model, where the features like the slope change or 
the low-lying continuum can be observed "generically" , 
without appealing to the existence of any additional in- 
teractions. However, one may in principle hope that sim- 
ilar effects could be observed in other one-dimensional 



G 



spin systems as w< 
mental indications^ 
present in CuHpCl. 



11. There are, for example, experi- 
I that such additional terms are really 



In recent experiment! on the S = 1 Haldane chain 
compound NDMAP it was observed, that one of the 
ESR branches was just continuing into the critical re- 
gion, without any noticeable features at H = H c . It 
is worth noting that this feature closely resembles our 
conclusions for the \a) — > |6) processes in the planar 
chain: if one denotes the (field-dependent) magnon gaps 
for NDMAP as A± and A , then at H < H c the above- 
mentioned ESR branch corresponds to the thermally ex- 
cited transition within the magnon triplet, with the fre- 
quency uj = Aq — A + . At H > H c there is a finite density 
of A + -magnons in the ground state, and, if we exploit the 
analogy with our picture of \a) — > \b) processes for the 
planar chain, one can expect the presence of a quasipar- 
ticle peak at exactly the same frequency, in agreement 
with the experiment. 

Finally, some words of caution are here in order. In 
the present paper, we have studied only the ground state 
transitions, in other words the ESR response at zero 
temperature, and only in purely one-dimensional (ID) 
model. When interpreting the experimental data, one 
should have in mind that, because of the gapless nature 
of the ground state in the critical phase of the ID sys- 
tem, the effects of temperature become important, as well 
as those of weak 3D coupling; particularly, the system 
should be 3D ordered under certain critical temperature. 
Thus, the results presented here should only be taken 
as a guide displaying features of the purely ID behav- 
ior. Further work is required to analyze the possibility of 
observation of those phenomena in realistic systems. 
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h c h s h 

FIG. 1. Schematic picture of ground state transitions in the ESR spectrum of the planar S = 1 chain as functions of the field 
in the critical region. Dashed line corresponds to the quasiparticle peak determined by a — ► b processes which are forbidden in 
absence of the in-plane anisotropy E. Filled areas show the continua, and solid lines within the continua indicate the position 
of the singularities. Arrows on the lines indicate the direction of the intensity increase when varying the magnetic field. 
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